Imaging of quantum Hall states in ultracold atomic gases 
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We examine off-resonant light scattering from ultracold atoms in the quantum Hall regime. When 
the light scattering is spin dependent, we show that images formed in the far field can be used to 
distinguish states of the system. The spatial dependence of the far-field images is determined by 
the two-particle spin-correlation functions, which the images are related to by a transformation. 
Quasiholes in the system appear in images of the density formed by collecting the scattered light 
with a microscope, where the quasihole statistics are revealed by the reduction in density at the 
quasihole position. 
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Electrons that are confined to two dimensions in a mag- 
netic field and cooled to within a few degrees of absolute 
zero behave in ways that have fascinated physicists for 
decades. In 1980 von Klitzing et al. performed a ground- 
breaking experiment with a cold electron system and ob- 
served plateaus in the Hall resistance as the magnetic 
field was varied Remarkably, the Hall resistance was 
exactly quantized at these plateaus to a value of h/(ne 2 ), 
where n is a positive integer. This integer quantum Hall 
effect is understood in terms of noninteracting electrons 
and is due to complete filling of the single-particle elec- 
tron orbitals known as Landau levels. 

This discovery was soon followed by the observation 
of a similar effect, where the Hall resistance was exactly 
quantized as h/(ue 2 ), with v a fraction This frac- 
tional quantum Hall effect occurs when the Landau levels 
are partially filled, with filling factor v, and the interac- 
tions between the electrons drive the gas into a strongly 
correlated state. This effect cannot be understood using 
perturbation theory and much of the current understand- 
ing is due to proposal of trial wave functions, such as the 
celebrated Laughlin wavefunction Q , given by 
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The integer and fractional quantum Hall effects were 
originally observed in electron systems where the elec- 
tron spin was polarized as a result of the Zeeman energy 
associated with the magnetic field. However, Halperin 
found that in some cases the Zeeman energy would be 
smaller than the other energy scales involved and that 
low-energy states with reversed spins could exist at some 
filling factors Halperin proposed a spin-singlet state 
with half the electrons in each spin state as a poten- 
tial ground state of the system. Halperin's intuition was 
confirmed by numerical studies of small electron systems 
that showed at filling factor v = 2/5 for vanishing Zee- 
man energy the exact ground state is unpolarized and 
has good overlap with the Halperin state [f| Q . 

Currently experimental studies of quantum Hall states 
use semiconductor materials that have their properties 
denned and limited by their composition. This makes it 



difficult to study the quantum Hall effect as system pa- 
rameters are varied. A more flexible realization of these 
effects may be possible in ultracold, neutral atomic gases 
@1 ■ Synthetic magnetic fields can be induced for these 
neutral gases by rotation Q or by optical means [l^, El , 
potentially providing access to the quantum Hall regime. 
It is hoped that the purity of systems of ultracold atoms 
and the ability to dynamically control system parameters 
will enable new insights into the nature of the quantum 
Hall effect and allow novel states to be achieved, such as 
those with non-Abelian excitations that may be used for 
topological quantum computing [12j. 

The usefulness of simulating the quantum Hall effect in 
ultracold gases will be limited by our ability to gain infor- 
mation about states of the system, and effective means 
of probing the system must be developed. A number 
of methods have been suggested. Characteristics of the 
wave function may be revealed by time-of-flight imaging 
[L3l4l5l ] or density profile measurement [l6| . Excitation of 
edge states reveals information about the systems topo- 
logical order [T3, EH an d Bragg scattering has been pro- 
posed to further probe the excitations of the system 
|20j. Density measurement in Fermi gases could be used 
as a probe of the effective conductivity [2lH23j and the 
statistics of quasiparticlcs might be probed in a Ramsey- 
type experiment 2J|. Here we examine what imaging 



based on light scattering reveals about an ultracold gas 
in the quantum Hall regime while the atoms remain in 
situ. 



I. LIGHT SCATTERING FROM 
MULTI-COMPONENT STATES 

We consider a specific example of light scattering from 
a Bose gas, where the N atoms are divided equally into 
two rriF spin states of a particular hyperfine level, which 
we refer to as spin up and spin down. In neutral atomic 
gases, synthetic magnetic fields can be created that only 
couple to the motional state of the atom and do not cou- 
ple to the spin of the atom, unlike an actual magnetic 
field. This extra experimental degree of freedom could 
allow the quantum Hall regime to be realized with zero 
Zeeman effect. For a two component Bose gas, with an 
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atomic interaction that does not depend on the internal 
state, a candidate grou nd state is the bosonic version of 
the Halperin state [25|, given by = S[ip], where 



N/2 



N 



k<l N/2<r<s 
N 

xn^-%K E » w2/j (2) 

and S symmetrizes over the particle labels. This state 
is expressed in terms of Zj = (xj + iyj)/£, where Xj = 
(xj,yj) are the coordinates of the atoms, and fjj 4-j are 
the spin-up and spin-down spinors. The characteristic 
length scale I = \J H/(eB) is determined by the effective 
magnetic field B generated by the optical field or rotation 
and is of the order of a few microns for typical param- 
eters [26]. This bosonic Halperin state has associated 
fractional filling factor v\ = 2/3. 

When the atomic interaction becomes spin dependent, 
such that the interaction between like spins is much 
stronger than the interaction between opposite spins, a 
paired state can exist [251 ] . Atoms of opposite spin pair 
together to form spin-triplets, resulting in a Pfaffian state 

*/• PI' r titj+Uj ^l flK -zOe-EJ-lV*, (3) 
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where the Pfaffian is defined for an N x N antisymmetric 
matrix M by [27j 
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This Pfaffian state is closely related to the ground state 
for the BCS theory of electrons with spin-triplet p-wave 
pairing and represents a superconducting state of effec- 
tive fermions [251 ] . 

Scattering of light has been suggested theoretically and 
used experimentally in a number of contexts to probe 
the spin distribution of ultracold atoms (28l - [39| . Here 
we wish to determine if light scattering is a useful probe 
of the many-body states described above. We consider 
experiments done with alkali-metal atoms where the de- 
tuning from the Sin P1/21P3/2 transitions is larger 
than the natural line widths and also larger than the 
Rabi frequency of any particular radiation mode. In this 
case the evolution of the atomic excited states can be adi- 
abatically eliminated and in the rotating wave and dipole 
approximations the light-matter interaction is described 
by 

H= /«L 2 W(E-(r)i lim! E+(r))f mi (r), 

mi ,rri2 
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where a mijm2 is the polarizability tensor (4fj, [41]. Here 
the atomic field operators ^^(r) and & m { r ) create and 
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FIG. 1. (Color online) Light propagating in the z direction is 
incident on the two-dimensional gas confined to the xy plane. 
Light scatters and is collected by an imaging system on the z 
axis. 



destroy atoms in the ground hypcrfinc state \Fm), and 
E + (r) and E~(r) are the slowly varying positive- and 
negative-frequency components of the electric field. 

The polarizability tensor can be decomposed into three 
parts ]40| , which correspond to a spin-independent inter- 
action that does not play a role in our imaging follow- 
ing the filtering procedure discussed below and two spin- 
dependent interactions. The second spin-dependent term 
goes to zero as the detuning is increased and is typically 
at least an order of magnitude smaller than the first and 
we may neglect it. We are then left with the following 
effective interaction [32, uSal 
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where p(r) = £ mi ,m 2 *L 2 (r)(Fm 2 \F/h\F mi )V mi (r) 
depends on the spin density of the atoms and a (A) 
is the detuning dependent coupling (see Hammerer et 
al. [40(). We assume in our example calculation that the 
two atomic spin states have tuf = ±1 and that the spin 
quantization axis is z. 

We consider a two-dimensional atomic gas in the xy 
plane that is illuminated by coherent light propagating 
in the z direction, wave vector kgz, with linear polariza- 
tion in the y direction as shown in Fig. Q] The inter- 
action described by Eq. (|6|) then leads to scattering of 
light to wave vectors k = ^(sin^ cos <j>, cos#, sin# sin0), 
where 9 is the polar angle with respect to the y axis 
and 4> is the azimuthal angle in the xz plane. The off- 
resonant nature of the scattering results in frequency 
shifts for the light of order of the atomic recoil frequency, 
which is many orders of magnitude less than the car- 
rier frequency, and to a very good approximation we 
have k — ko. The wave vector of the scattered light 
it then completely determined by k, the projection of 
k into the xy plane. For light scattered predominately 
in the z direction the interaction results in a rotation of 
the polarization, where we describe the scattered light 
in terms of polarization vectors £i(k) = (sin (j), 0, — cos <f>) 
and £2 («0 = (— cos^ cos <fi, sin 6?, — cos# sin</>). Placing a 
polarization filter with extinction axis in the y direction 
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FIG. 2. (Color online) Far-field images of (a) the bosonic 
Halperin state, (b) the Pfaffian state, and (c) the compar- 
ison Laughlin state, (d) Slices through the far-field images 
at K y = for the bosonic Halperin state (solid), the Pfaffian 
state (dashed), and the comparison Laughlin state (dotted). 
Images are generated from 100000 Monte Carlo configurations 
of 100 particles. 



on the z axis will then remove the input light from the 
signal and will also filter the majority of light scattered 
into the mode with polarization s2(k), as well as light 
scattered by the spin independent interaction [32|. 

Image formation then depends on the rate of scattering 
to k with polarization ei(/c). To first order in the cou- 
pling constant [42j . this is proportional to I(k) — ((f\<p) 
where 
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for the many-body state and the integration has been 
reduced to two dimensions by taking into account the 
tight confinement in the z direction. The first-order 
treatment of light scattering is expected to reproduce 
the same images as an experiment, provided the ex- 
perimental images are averaged over a number of runs 
|43| . For scattered light propagating predominately in 
the z direction, ci(k) X y = zsin0 + xcos0 ~ z and 
we can neglect the small corrections resulting from the 
x component |38[. The dominant features of an image 
then result from the z component of the spin operator, 

One of the simplest methods of imaging the sys- 
tem is to let the scattered light propagate into the 
far field, where the distance from the sample is much 



greater than the wavelength of the imaging light. In 
this region the spatial dependence of the light in- 
tensity is simply determined by the Fourier compo- 
nents of the scattered light, that is /(ax) oc 
with scaling constant a that depends on the distance 
from the sample to the detector and determines the 
spatial size of the image. For a many-body state 
I*} = JT I rfxi---dxAr*(xi,--- ,xjv)|xi, • • • ,xjv) where 
the first N/2 coordinates refer to spin up atoms and the 
second to spin down, we have 



= J dx!---dx N \^(xi,--- ,xjv) 
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The far- field images for the states from Eqs. ©-© can 
then be calculated by using the Metropolis Monte Carlo 
method to evaluate the integral over N coordinates (see, 
e.g., Ref 0). 

In Figs.[2ja) and[2jb) we show the far field image of the 
bosonic Halperin state from Eq. ([2]) and the Pfaffian state 
from Eq. ([3]). For comparison in Fig. [5Jc) we show the 
image of a state with N/2 atoms of each spin but with 
the spatial wave function given by the Laughlin wave 
function in Eq. ([1} with n = 2. The images are generated 
by averaging over 100000 Monte Carlo configurations of 
100 particles. The far-field images allow us to distinguish 
the three different states, where all the images go to zero 
at k = but the nature of this approach varies for each 
state as shown in Fig. ^d). The difference is due to 
the different spatial correlations between opposite spins 
in each wave function. In the Halperin state the atoms 
of opposite spin are allowed to exist closer to each other 
than particles of the same spin, which results from the 
lower power of the (zj — zj) terms for atoms of opposite 
spin in the wave function. This leads to long wavelength 
oscillations in the far field that create the large dip in the 
image of this state. For the Laughlin state the particles 
of opposite spin interact in the same way as particles of 
the same spin and these long-wavelength oscillations do 
not occur. In the Pfaffian state the (zj — zi) term does not 
enter the wave function for atoms in a pair, and paired 
atoms come much closer to one another than other atoms 
leading to an intermediate dip at n = 0. 

In fact, the far-field image of the gas is related to its 
two-particle spin-correlation function by Fourier trans- 
form. We have 
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are the two-particle correlation functions, which give the 
probability of finding a particle in state a separated by 
x from a particle in state (3. In Fig. [HJa) we plot a slice 
through the transform of the far-field images given by 
Eq. ©. In Figs. EIb)-[3](d) we show the corresponding 
correlation functions for the Halperin state, the Pfaffian 
state, and our comparison Laughlin state. The combi- 
nation of the correlation functions given in Eq. © leads 
largely to cancellation in the Laughlin case as all the cor- 
relation functions are equal. In the Halperin case g^(~x) 
exhibits extra lobes at x ~ 2.5£ that do not cancel in the 
sum and show up in the transformed image, similarly in 
the Pfafhan case (744- (x) is non-zero close to x = which 
does not cancel with the other terms. 

In Figs. [5] and [3] we have averaged over 100000 Monte 
Carlo configurations to estimate the expectation value 
in Eq. (|8]l. In a real experiment an image from a sin- 
gle experimental run will not give the expectation value, 
the combination of scattering and detection of photons 
will disturb the state leading to relative localization of 
the atoms (43j . and the image will be determined by the 
equivalent of a single Monte Carlo configuration. This 
does not mean 100000 experimental runs are required to 
distinguish the states; in fact, images formed by aver- 
aging over 10 configurations, although noisy, show the 
characteristic differences between the states allowing the 
states to be distinguished in a realistic number of experi- 
mental runs. Furthermore, because of the radial symme- 
try of the many-body state finding the radial intensity 
of the image by averaging over all angles can improve 
the ability to distinguish the states for a small number 
of experimental runs. 

Another experimental issue is the effect of the trapping 
potential used to trap the ultracold gas on the many- 
body state. For optically induced gauge potentials the 
trap potential and optically induced scalar potential can 
be arranged to cancel [HI • For rotating gases in the high 
rotation limit the trapping potential is a crucial factor 
in realizing the effective Hamiltonian of a magnetic field, 
and at the critical rotation frequency the trapping poten- 
tial plays no additional role. In both these experiment 
arrangements the many-body states of th e sy stem are ex- 
pected to be the ones considered above [25|; however, if 
an excess trapping potential remains in addition to the 
magnetic-field Hamiltonian, the state can form domains 
of quantum Hall states with different filling factors that 
would lead to an altered far- field image [16|. We do not 
consider this case here. 



II. IMAGES OF QUASIHOLES 

It is also interesting to consider what imaging can re- 
veal about quasiparticles in these quantum Hall systems. 




FIG. 3. (a) Slice along y = of the transform given in Eq. ([9]) 
of the far-field images for the Halperin state (solid) , the Pfaf- 
fian state (dashed) and the comparison Laughlin state (dot- 
ted). [(b)-(d)] Two particle correlation functions for the (b) 
Halperin state, (c) the Pfaffian state, and (d) the Laughlin 

state, with g++(x,y = 0) = g (x,y = 0) (solid lines) and 

g+-(x,y = 0) = g- + (x,y — 0) (dashed lines). Functions 
were calculated using 100000 Monte Carlo configurations of 
100 particles. 



For multicomponent systems, quasiholes can be created 
by a focused laser which repels atoms at a particular posi- 
tion, say, rj. The laser can be arranged to primarily affect 
only one of the spin components or, alternatively, mul- 
tiple components [25j . For the two-component bosonic 
Halperin state where the the laser affects only the spin- 
up atoms the state is then given by 



4L 



S 



'jV/2 

IK* 

i=l 



(11) 



5 



(a) (b) 




FIG. 4. Images of the two-component Halperin state with 
a spin-up quasihole at (x/i,y/£) — (5,0). (a) Total density, 
(b) spin-up density, and (c) spin-down density, (d) Slice with 
y/t = through images of the (solid) total density, spin- up 
density, (dotted) and spin-down density (dashed). Images 
were calculated using 100000 Monte Carlo configurations of 
100 particles for an imaging system with numerical aperture 
sin (9 = 0.5X/1. Images are normalized to the total number of 
particles. 



can be used to reveal the density of the state [38j. In 
this case, light scattered by the sample is collected by a 
microscope with numerical aperture sm8 and forms an 
image of the density or spin density. We can choose the 
probe light so that it scatters identically from either spin 
state or so that it scatters from only one state, giving a 
total density image or an image of the density of either 
spin state. 

Images of a bosonic Halperin state with a spin-up 
quasihole at (x/£,y/£) — (5,0) are shown in Fig. H) We 
can see a dip in the total density and a larger than nor- 
mal density of spin-down atoms at the quasihole position. 
In fact, the quasihole statistics are revealed from these 
images by integrating the image intensity over the region 
of excess or reduction to find the volume added to or re- 
moved from the normal density. We find that the dip in 
the total density has a volume of 1/3 of an atom, and the 
peak in the spin-down atom density also has a volume of 
1/3 of an atom. The dip in the spin-up density has a 
volume of 2/3 of an atom. This property is also seen 
in the images of similar states with more than two spin 
components. For the n-component generalization of the 
Halperin state, the quasiparticles have p/(l + n) statis- 
tics, where p is the number of components the quasihole 
laser affects [25[ . Then, for example, in images of a three- 
component gas where two components arc affected by the 
laser, the total density image has a dip of volume 1/2 of 
an atom at 77, while the unaffected spin component has a 
bump of volume 1/2 of an atom at 77. In general, the to- 
tal density and the unaffected components will have less 
and excess, respectively, of p/(n + 1) of an atom at 77 and 
the affected components will have 1 — p/(n + l) less of an 
atom at 77. 



III. SUMMARY 



In effect, the spin- up particles see an extra spin-down 
particle at 77, leading to a repulsion from 77 and a lower 
number of spin-up particles in the region close to 77. The 
spin-down particles, on the other hand, are not affected 
by the laser and a higher than normal density of spin- 
down particles will exist at 77. This type of quasihole is 
an anyon with 1/3 statistics. 

Far-field images of states with a small number of quasi- 
holes do not differ significantly from the corresponding 
image of the state without quasiholes. This is because 
adding a quasihole causes only a local defect to the state 
and hence does not change the particle-particle correla- 
tion function significantly. Instead, microscopic imaging 



The above examples demonstrate the potential that 
imaging of scattered light has as a tool for investigating 
the properties of ultracold atomic gases in the quantum 
Hall regime. Imaging in the far field provides a probe of 
the density-density correlation of the atomic matter and 
when the light scattering is spin dependent it provides 
a probe of the spin-density correlations. Different states 
then have different optical signatures that relate directly 
to the spin-correlation functions of the state, enabling 
the identification of particular states in experiment. Fur- 
thermore, imaging the density using a microscope can 
reveal quasiholes, where the reduction in the density at 
the quasihole position gives the quasihole statistics. 
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